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THE DYNAMIC MODELLING OF A SLOTTED 
TEST SECTION 


1. Introduction . 

In designing the control system for the National Transonic Facility, 
it was necessary to develop a mathematical model of the wind tunnel dynamics. 

Due to the limitations of present day modelling, it was necessary to restrict 
the modelling techniques to the use of one dimensional unsteady flow. With this 
limitation it was still imperative that the dynamic characteristics of a slotted 
test section be incorporated into the model. This report describes the model 
developed to accomplish these objectives. 

2 . Objectives . 

The following objectives were set for the test section-plenum model to be 
developed. 

a. The model should serve as a means for determining the dynamic 
performance of the NTF and the subsequent design of optimal control configurations. 

b. The model should be of the one dimensional flow type so as to 
facilitate computer simluation. 

c. The model should allow for the simulation of the effects of test 
section blockage, boundary layer losses, slot flow required for supersonic flow 
generation, reentry flap position, and auxiliary plenum suction. 

d. The model should be capable of validation by comparison to previous 
tests made on the eight foot tunnel at Langley. 

e. The model should be capable of being updated as operating data from 
the NTF becomes available. 



f. The model should lead heuristically to simpler models so 
that the principles of operation of this part of the tunnel are more clearly- 
understood. 

The model presented in this report meets most of the objectives. 

The shortcomings of the current representation are: 

a. There is no representation of reentry flap position. 

b. The boundary layer coefficients Cxm are not known to a 
sufficient degree of precision. 

c. The test section divergence is not represented. 

d. There is no simplified representation available so as to 
better understand the test section dynamics. 

e. The actual value of Ar is not known. Ar is an equivalent area 
and may not be equal to the geometric area. Ar may even vary with mach 
number. 
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f. The assumption of a constant pV loss coefficient for the 
diffuser may not be a good one. 

3 . Basic Configuration . 

The basic configuration of the slotted test section and plenum is 
illustrated in figure 1. The following equations are used: 

3.1 The flow entering the test section is given by: 


Qthrt = 


^ 144 •Atest*Pstag»Mach 

(1 + .Z.Mach^) ^ 


144* Atest-Pstag 

( 1 . 2 )^ 


1.4. g 

R«Tstag 


if Mach < 1 


( 1 ) 


1.4. g 


R*Tstag 


if Mach > 1 


* see Appendix B 
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where: 
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Atest is the test section area, Atest = 66.77 ft 

R is the universal gas constant, R = 55.1 ft/®R 

g is the gravity constant, g = 32.2 ft/sec^ 

Pstag is the stagnation pressure (Ib/in ) 

Tstag is the stagnation temperature (°R) 

Mach is the test section mach number 

Pstag and Tstag are determined by the tunnel simulation. 
3.2 The mach number is given by: 


Mach = 



2/7 

fPstag 1 


LPstest J 


5 


( 2 ) 


where: 
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Pstest is the static test section pressure (Ib/in ) 

3.3 The flow (Ibs/sec) entering the plenum from the test section via the 
slots is given by: 


Qslot = QthrfCxm + (Qthrt-Qtest) + Qdm (3) 

where: 

Qb = Qthrt»Cxm is the component due to boundary layer. 

Cxm is a function of mach number (see section 3). 

Qt = (Qthrt-Qtest) is the component due to supersonic flow 
generation. 

Qdm = Qthrt*Qdmx is the component due to model blockage. 

Qdmx is an input variable. 



3.4 The flow exiting the test section is given by: 


Qexit = Qthrt - Qslot 

3.5 Qtest is the test section flow and is given by: 

144* Atest«Pstag*Mach / 1.4*g 

Qtest = 2 — T '\/ (^) 

(1 + .2‘Mach ) V R'Tstag 

3.6 The flow leaving the plenum and flowing into the diffuser via 
the reentry flaps is given by: 

Qflap = Qr - Qexit (5) 


where: 

Qr is the total flow into the diffuser. 


3.7 The total weight of gas in the plenum is given by: 


Wpl = 


(Qslot - Qflap) dt 


3.8 The static pressure in the plenum is given by: 


Ppl = 


Wpl-R*Tpl 

144-Vpl 


( 6 ) 


(7) 


where: 

Tpl is the plenum static temperature. It is assumed 

that the total and static temperature in the plenum 
are equal to each other as well as to the test section 
total temperature. 
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Vpl is the plenum volume, Vpl = 36,000 ft 
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3.9 The test section static pressure is assumed equal to the plenum 
static pressure. 



Pstest = Ppl 

(8) 

3.10 

The static temperature in the test section is given by: 



Ttest 

Tstest = — 

(1 + .2 Macli ) 

(9) 


Ttest is equal to Tstag less any thermal losses. 


3.11 

The total temperature in the plenum is given by: 



Tpl = Epl/(Wpl*Cv) 

(10) 


• 

Epl = Cp(Qslot*Ttest - Qflap«Tpl) 

(11) 

3.12 

The velocity in the test section is given by: 
Utest = Mach V' 1.4*g*R*Tstest 

(12) 


3.13 The ejector action of the diffuser is represented by: 

144*Ar*Pstest + Qexif Utest/g = 144*Ar*Pr (1 + 1.4‘Mr^) (13) 

where: 

The balance is taken across an equivalent area Ar. The exact 

value of this area is not known. A nominal value of 
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“ 1.16 Atest = 77.45 ft was chosen. The value of Ar was varied 
during the simulation so as to obtain a measure of the sensitivity of 
the system to the value of Ar, This is the familiar momentum equation 
applied across a line. The equation balances the sum of the potential and 
and kinetic forces on either side of a line. On the left face of Ar, the 
potential force is 144*Ar •Pstest . The kinetic force is given by 



Qexif Utest/g, since we may assume that the kinetic energy associated with 
the flow Qflap is negligible. The energy on the right face of Ar is 
obtained by using the identity. 
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144*Ar*Pr*g 1.4 Mr = Qr*Ur (14) 


3.14 The remainder of the tunnel, up to the fan, is represented by a 
sequence of momentum, energy, and continuity equations. The first 
equation encountered is: 


^9 = 


144*A. 


■g (Pr - - Xq^Q^Ug) + QrUr - 


/XL, 


(15) 


where: 

The equation is written across a section of the diffuser of 

2 

length XLg, XL^ = 30ft, and an equivalent area of Ag, Ag = 90ft . 

The state of the gas entering the section is (Pr,Qr,Ur). The state of 
the gas leaving is Qg assumed equal. This 

assumption is valid since it is equivalent to saying that the flow 
entering a volume of short length is equal to the average flow in that 
volume. 
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Xq^. is the pV loss coefficient in this section. 

Xqg = .1111 10 ^ sec^/(in^f t) . 

4 . Solution of Equations 

The equations of 3.1 to 3.14 are solved simultaneously in an 
iterative manner. Care must be taken in the calculations so as to remain 


on the proper branch of the solution. This will assure that choking and 
the shock equation are properly represented. The procedure used is out- 
lined in Appendix A. In general, the procedure proceeds as follows: 
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Let 


Fm = 144*Pstest* Ar + Qexit*Utest/g 


(16) 


and 


d = Qr^R Ttr /(1.4‘Fm^‘g) 


(17) 


whsrG Ttr is the total temperature on the right face of Ar. If 
there is no loss of heat in the plenum Ttr = Ttest. Mr is given by the 
quadratic expression. 

(5-49 d)*Mr^ + (25 - 70d)*Mr^ - 25 d = 0. (18) 

The remaining variables are given by: 


Pr = Fm/(144*Ar(l + 1.4Mr^)) 

Tr = Ttr/(1 + . 2Mr^) 

Ur = 49.84*Mr ~\fTc 

If Mr > 1, choking occurs and we set Mr equal to 1. 
Qr is then given by: 


Qr = Fm y^5(1.4)g/(24-R*Ttr) 


Qg is set equal to Qr, and Qg is constrained so that Qg < 0. 

5. The Distributions of Losses Around the Tunnel . 

For simulation purposes, the tunnel has been divided into 
15 stations. Where: 

Station 1 is at the fan outlet 

8 is at the stagnation chamber 
10 is at the entrance to the high speed diffuser 
15 is at the fan inlet. 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 


From the report by D. M. Rao (11), it may be inferred that when 
operating at 14.7 psia, 500°R the tunnel drops will be distributed as 
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follows 


station 

1 to 4 

13.91% 

station 

4 to 6 

3.21% 

station 

6 to 8 

11.37% 

station 

8 to 10 

60.86% 

station 

10 to 12 

5.46% 

station 

12 to 14 

5.19% 


From the report by B, Gloss (2) the fan ratio as a function 
of mach number is given as: 


Mach No. 

Fan Ratio 

.2 

1.0120 

.3 

1.0180 

.4 

1.0245 

.5 

1.0330 

.6 

1.0430 

.7 

1.0550 

.8 

1.0700 

.9 

1.0880 

1.0 

1.1120 

1.1 

1.1420 

1.2 

1.1850 


5.1 Using the above data the losses around the tunnel may be established 
if we further assume: 
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a. The losses are pV losses, with constant loss coefficients. 

b. That the slot flow due to boundary layer at Mach 1 is 
3% of the flow at the throat. 

c. That the slot flow is a function of mach number. 

5.2 The coefficients i = 3,5,7,9,11,13 are defined by: 


DP. 

X 


Xq . • Q . • U . 

^x X X 


(23) 
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where: 


is the pressure loss in the ith section, 
is the loss coefficient in the ith section, 
is the flow in the ith section, 
is the velocity in the ith section. 


The coefficients Cxm^ are defined by: 


Qb = Cxm^*Qthrt 


where : 


The subscript i denotes the mach number at which the 
computation is made. 

Qb is the slot flow due to boundary layer. 

Cxm^ is the slot flow coefficient. 

5.3 The coefficients Xq and Cxm are determined as follows: 

a. Establish the tunnel simulation at 14.7 psi, 500°R, 

Mach 1.0. Fix Cxm at .03 and vary all of the ^0^ in an iterative fashion 
the fan ratio is 1.112 and the losses are distributed as indicated 
in 5.0. 


b. Establish the mach number at .1 to 1.2 in steps of .1 and 
for each mach number vary Cxm^ until the proper fan ratio is obtained. 
If this is done, the following values are obtained: 


Xq3 = 

.8055 

10 ^ 

Xqg = 

.3491 

10-' 

X 

II 

.1934 

10-^ 

Xqg = 

.1111 

10-6 



-7 

Xqii = 

.1791 

10 

Xqi3 = 

.2013 

10-' 
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Mach No. 


Cxm 


.1 

.4501 

.2 

.1612 

.3 

.0807 

.4 

.0445 

.5 

.0300 

.6 

.0232 

.7 

.0206 

.8 

.0216 

.9 

.0243 

1.0 

.0300 

1.1 

.0299 

1.2 

.0216 


5.4 The data given in 5.3 is suspect due to the large values of Cxm 
at the low mach numbers. This is more than likely due to the inaccuracies 
of the assumed fan ratios. In order to correct for this error we introduce 
the energy ratio. 


Qthrt Utest^ 

CE = (25) 

2(550)g HP 


where HP is the horsepower needed to drive the tunnel. 

The coefficients Cxm are then redetermined. Below mach numbers of 
.6 the fan ratios as given in 3.0 are ignored. Instead the values of Cxm 
are chosen so that the value of CE is held constant and equal to that at 
Mach = .6. This results in coefficients which appear to be more realistic. 
It is these values that are to be used to represent the test section. The 
results of this calculation is tabulated on the following page. 
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Mach No. 

Cxm 

CE 

Cr 

Crl 

.1 

.0137 

4.70 

1.0012 

1.0080 

.2 

.0144 

4.70 

1.0048 

1.0120 

.3 

.0155 

4.70 

1.0107 

1.0180 

.4 

.0169 

4.70 

1.0189 

1.0245 

.5 

.0210 

4.70 

1.0299 

1.0330 

.6 

.0232 

4.70 

1.0430 

1.0430 

.7 

.0206 

4.73 

1.0550 

1.0550 

.8 

.0216 

4.96 

1.0700 

1.0700 

.9 

.0243 

4.89 

1.0880 

1.0880 

1.0 

.0300 

4.65 

1.1120 

1.1120 

1.1 

.0299 

4.39 

1.1420 

1.1420 

1.2 

.0215 

4.01 

1.1850 

1.1850 


Where Cr is the resulting fan ratio and Crl is the fan ratio given in 5.0. 

6.0 The time response of the test section . 

6.1 It was suspected from the system equations that the dynamics of the 
test section can be represented by a simple first order system. To 
verify this, the time response of mach number due to a 2 degree step change 
in guide vane angle was determined by the NTF dynamic model. The quantity 



where Mach, Mf, and Mi are the current, final, and initial values of mach 
number, was plotted as a function of time on figure 5. The quantification of 
the points is due to the fact that Mach, Mf, and Mi were only observed to 3 
significant digits. 

Since a straight line fits through these points it is verified that 
the system is of the first order. 

The X intercept of .35 corresponds to a dead time of .35 seconds. 

This agrees with the acoustic propagation time of 

343 ft/(863 ft/sec) = .397 sec (27) 

It may be shown that the slope of the line is the reciprocal of the 
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time constant. This yields a time constant of 1/.411 = 2.43 sec. If the 
time constant is determined in the usual fashion, by defining the time 
constant as the time required to reach .632 of the final change, a value of 
2.8 seconds is obtained. Note that 2.43 + .35 = 2.78. 

From this we conclude that one can model the test section response 

by: 

, -Td*S 
k*e 

AM = AG (28) 

(Tc*S + 1) 

where: 

J^G is the change in guide vane angle. 

M is the corresponding change in Mach number. 

Tc is the time constant. 

Td is the dead time, 

k is the gain. 

S is the Laplace transform operator. 

In designing the mach number control system, the above parameters 
prove useful. In the next section we examine how these parameters vary with 
the operating point. It is seen that they are a function of pressure, mach 
number, and Ar. In addition, a nonlinear affect is observed in that the time 
constant is dependent on the step size. Further, for small values of Ar 
choking occurs at supersonic operation. This choking occurs at the 
boundary of Ar and causes velocity limiting which in turn slows down the 
response. 

6.2 The variation of gain and time constant. 

On figures 2a and 2b the gain and time constant of the test section 
is plotted as a function of stagnation pressure and mach number. This 
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data was obtained by runs made on the NTF simulation. With the pressure 
and temperature control active a 2 degree step change in guide vane angle 
was made. The time constant was determined as the time required for 
the mach number to make .632 of the total change. Since we have already 
shown that the system is indeed a first order, this calculation is valid. 
However, the time constants shown are the true time constant plus the 
system dead time. Since the distance between the test section and the 
fan is approximately 343 ft, the dead time varies between .307 and .487 
seconds for the temperature range of 500 to 200'^R. To obtain this data, 
runs were made at pressures of 30, 60, 90, and 120 psi; temperatures of 
200, 300, 400, and 500°R; and mach numbers of .2, .3, .4, .5, .6, .7, .8, 

.9, 1.0, and 1.1. There was no change in either gain or time constant with 
temperature. Since it is known that the dead time does vary with tempera- 
ture, it must be true that the actual time constant also varies with 
temperature in such a manner as to maintain the apparent time constant (as 
plotted in fig. 2a) fixed. The gain plotted on fig. 2b was obtained by 
taking the ratio of mach number change and guide angle change and multiplying 
by 100. 

6.3 The time constant variation due to Ar. 

It was observed that the value of Ar had a strong effect on the time 
constant. The curves of figure 4 sho^^^ the variation of time constant with 
Ar. It is noted that for low values of mach number that the value of Ar 
has little if any effect on the time constant. The time constant for 
transonic operation is strongly dependent on Ar. These curves clearly show 
the need for obtaining better information as to the value of Ar. 
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6.4 The time constant due to blockage. 

On figure 3, the time constant due to a 2% blockage is plotted. 
Comparing this with figure 2a, it is seen that the time constant for 
blockage and guide vane deflection is almost the same. That this is so 
is clearly seen by examining figures 5 and 6. From these we conclude 
that the response due to blockage differs from that due to guide vane 
deflection only in that blockage does not have any dead time associated 
with it. Note that figure 5 has an x intercept of .35 seconds, while the 
line of figure 6 goes through the origin. 

6.5 Comparison with the eight foot tests. 

In order to assist in the validation of the NTF simluation transient 
tests were made on the eight foot tunnel at Langley. These tests were made 
in order to determine the response of the plenum-test, section blockage. 

The findings of these tests are documented in ref. (10). The table below 
which appears on page 11 of the above report summarizes the pertinent 
results. 


Test 

Run 

Nominal 
Mach No. 

SEL (#44) 

Average Time Constant, sec 
FM(IOIT) Eqs. (11) to (13) 

19 

1.05 

7.7 

9.2 

13.5 

20 

1.20 

3.7 

5.2 

3.5 

21 

.95 

7.1 

6.7 

4.5 

22 

.80 

5.3 

5.5 

4.6 


In order to compare the present NTF simulation with these tests, 
the NTF simulation was altered so that it would exhibit the properties of 
the eight foot tunnel. This was accomplished by changing volumes and lengths 
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in the tunnel proper. This change is described in ref. (8). The plenum 

q 

volume was increased from 36,000 to 72,000 ft . The value of Ar was left 
2 

at 77.45 ft , and a linear fan representation was utilized. The inbleed 
and outbleed valves were locked and a 2% step change in blockage was made 
in the test section. Four runs were made: run 2219, Tstag = 500, 

Pstag = 14.71, Mach = 1.05; run 2220, Tstag = 500, Pstag = 14.71, 

Mach =1.2; run 2221, Tstag = 500, Pstag = 14.71, Mach = .95; and run 2222, 
Tstag = 500, Pstag = 14.71, Mach = .8. These correspond to runs 19, 20, 

21, and 22 respectively. A typical response is shown on figure 8. 

An analysis of the output of figure 8, is made by the use of 
figure 7. From this figure we conclude: 

a. The response of mach number is of the first order. For run 2219, 

2 

with an Ar of 77.45 ft , the time constant is 3.95 see, there is no system 
dead time. 

b. The response of the test section pressure is of the first order. 
For run 2219, with an Ar of 77.45, the time constant is 3.98, there is no 
system dead time. 

c. The response of stagnation pressure, is of the third order. 

For run 2219, with an Ar of 77.45 ft“, there is a time constant of 3.69 
seconds, with an additional underdamped second order of approximately 

3 hertz. There is a system dead time of .35 seconds. 

The table below compares the eight foot tests with the simulation 
results . 
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Run 

8 ft tests 

simulation 
Ar =70 Ar = 72 

> 

II 

Ar = 77.45 

19 

7.7, 9.2, 13.5 

C 

6.3 

5.3 

4.0 

20 

3.7, 5.2, 3.5 

C 

C 

19.0 

2.9 

21 

7.1, 6.7, 4.5 

6.2 

5.0 

5.6 

3.9 

22 

5.3, 5.5, 4.6 

4.9 

4.6 

4.4 

4.7 


The listed values are measured time constants in seconds. A C 
indicates that the simulation failed to run due to choking at Ar. 

Examining the above, it may be concluded that the simulation 
concurs with the test. In order to obtain better correlation or detect 
differences between the simulation and the tests, it will be necessary to 
obtain more consistent test data and establish more precise values of Ar. 
7.0 Mach Number Control 

From the low band-pass of the mach number response, it is apparent 
that the guide vane actuator will be the lim.iting factor on the response 
of. the mach number control system. 

The process time constant is of the order of 3 seconds, while the 
actuator is a second order system of 6 radians with a damping factor of 
.8. Together these two yield a phase angle lag of 130° at 3 radians, which 
limits the response to 3 radians per second. This is in agreement with the 
analysis presented in (6) . 

The settling time, to within 1%, for a second order system with a 
natural frequency of 3 radians per second and a damping factor of .8, 
is obtained by solving 
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.01 


Vl - .82 

This yields a settling time of 2 seconds. 

This is not consistent with the system perforaance as exhibited by 
the NTF simulation. A typical simulation output is contained on 
figures 9a, 9b, and 9c. This run corresponds to a step change in mach 
number set point of .1, with the tunnel operating at 100 psi, 300®R, and 
mach = .8. In this run a settling time of 20 seconds was required. This is 
far different from that expected for a properly tuned system. The discrep- 
ancy is due to cross coupling, the presence of reset in the control, 
improper tuning and perhaps some nonlinear affect unaccounted for. This 
discrepancy is to be investigated in a later study, meanwhile the target 
values of 2 to 3 seconds for settling times of the mach number control 
system is to be maintained. 

8.0 Conclusions and Recommendations . 

From the results presented in this report, it may be inferred that 
the proposed test section math model is realistic. It is therefore reason- 
able to proceed with the dynamic investigation of the NTF utilizing the 
model established in this report. 

The following action items have arisen as a result of this study: 

1. Obtain better data as to the distribution of losses in the 
tunnel. 

2. Establish more precise values of fan ratios as a function of 
mach number. 

3. Obtain better data so as to establish the boundary layer losses. 

4. Obtain data as to the variable losses in the high speed 
diffuser. 
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5* Determine more precise values for the parameter Ar. 

Establish any variation in Ar with mach number. 

6. Design an experiment on an operating tunnel to facilitate 
further validation of the simulation. 

7. Proceed with the design of a suboptimal mach number 
control system. 
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TIME PSTEST PSTAG MACH 
0,00 7,5257 10,7100 1,0500T 
,50 7,0020 10,6975 l,0O0OT 

1.00 7,0788 10,6700 1,0500T 

1.50 7,5550 10,6600 1,0210 T 

2.00 7,6166 10,6505 1,0139 T 

2.50 7,6715 10,6015 1,0070 T 

3.00 7,7227 10'6300 1,0009 T 

5.50 7,7657 10,6260 ,9957 T 


0,00 7,8006 10,6188 ,9909 T M 

0,50 7,8385 10,6123 ,9868 T M 

5.00 7,6677 10,6060 ,9832 T M 

5.50 7,8903 10,6007 ,9800 T M 

6.00 7,9181 1015972 ,9772 T M 

6.50 7,9388 10,5959 ,9707 T M 

7.00 7,9572 10.5909 ,9725 T M 

7.50 7,9736 10,5881 ,9705 T M 

8.00 7,9885 10,5803 ,9687 T M 

8.50 8,0015 10,5816 ,9671 T M 

9.00 8,0122 10,5800 ,9658 T m 

9.50 8,0212 10,5787 ,9607 T H 

10.00 8,0296 10,5771 ,9637 T M 

10.50 8,0377 10,5757 ,9627 T M 

11.00 8,0009 10,5703 ,9619 T M 

11.50 8,0512 10,5733 ,9611 T m 

12.00 8,0565 10.5725 ,9605 T M 

12.50 8.0607 10,5717 ,9600 T M 

13.00 8.0600 10.5706 ,9595 T 

13.50 8,0680 10,5700 ,9591 T m 

10.00 8,0711 10,5697 ,9587 T m 

10.50 8.0738 10.5695 ,9580 T M 

15.00 8,0762 10,5691 ,9581 T M 

15.50 8,0786 10,5686 ,9578 T m 

16.00 8,0809 10.5678 ,9575 T M 

16.50 8,0825 10,5676 .9570 T m 

17.00 8,0836 10,5676 ,9572 T M 

17.50 8,0800 10,5675 ,9571 T M 

18.00 8,0856 10.5673 ,9570 T M 

18.50 8,0869 10,5670 ,9568 T M 

19.00 8,0881 10,5669 ,9567 T M 
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Appendix A. 


The solution of the diffuser-ejector equation 


1. The momentum balance across the area Ar is given by: 

144 Ar Pstest + Qexit Utest/g = 144 Ar Pr (1 + 1.4 Mr^) 


2. The momentum equation between Pr and P^^ is given by: 
Qg = [1^4-Ag-g (Pr - P^Q - Xq^-Qg-Ug) + Qr-Ur - 


3. To solve these, in terms of Pstest, Mach, and P^^, we write: 


Mr = 


Ur 


Qr 


4*R*g*Tr 


Ar • "Yl.4*R*g*Tr 


Qr. (5 + 7 Mr ) / Ttr. R 


Fm.5 (5 + Mr ) V 1.4*g 


where: 


Fm*144 .pstest *Ar 4- QexifUtest/g 


Letting 


Qr .R.Ttr 


d = 


1.4 *Fm . g 

the above equation may be written as; 


5.Mr^ (5 + Mr^) 
(5 + 7-Mr^)^ 


= d 


or 


(5 - 49 d) Mr^ + (25 - 70 d) Mr^ - 25 d = 0 


( 29 ) 

(30> 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 
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5. The above system is solved by letting: 


25 - 70 d 25 d 

b ; c 

2 (5 - 49 d) (5 - 49 d) 

which yields the quadratic 


( 37 ) 


Mr^ - 2b Mr^ + c = 0 
We consider three cases 
case 1 . d £ 5/49 

O 

Then b < 0 and c < 0 and since Mr ^0, the unique 
solution is 


(38) 


Mr^ = b + Vb^ - c ( 39 ) 

case 2 . 5/49 < d < 5/24 

2 

There are two real solutions for Mr . Moreover if x is 
2 

a solution for Mr , then y is the dual solution. With x and y 
satisfying: 


5 + y 5 + X 

X = ; and y = (40) 

7y - 1 7x - 1 

further if x ^ 1, then y £ 1, and visa versa. * 


In this case we chose the solution for Mr < 1 or 



For d ^ 5/24, the solution is degenerate. In this case we set 
Mr equal to 1. The value of Qr is obtained from: 


(41) 


* We note that X=l/7 is degenerate. From the shock tables ( NACA Report 1135 ) 
it is noted that the mach number upstream of a normal shock wave is assymptotic 
to the square root of 1/7 as the downstream mach number goes to infinity. 
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Qr 


Fm 


(42) 


f5*l*4*g 


24*R*Ttr 


4.0 The ejector equation is seen to be consistent with the normal 
shock equation, if Qr = Qexit. For this case we get due to the 
symmetry between Mach and Mr. 


5 Mr^ (5 + Mr^) 5 Mach^ (5 + Mach^) 

5+7 Mr^)^ (5+7 Mach^)^ 


(43) 


If Mach < 1, then d ^5/24 and either case one or case two holds. 
For case one Mr = Mach, since there is only one solution. For case tx^ro, 
again Mr = Mach, since the other solution 


„ 5 + Mach^ 

Mr = (44) 

7 Mach - 1 

is not physically realizable (Mr can not be greater than one) . 

If Mach ^ 1, then 5/49 < d _< 5/24, and case two holds. 

There are now two possible solutions 

Mr = Mach, and Mr = V(5 + Mach^)/(7 Mach^ - 1) 

which is the well-known shock equation. 
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APPENDIX B - The Reentry Area Ar 


Immediately downstream of the test section there is a step increase in tunnel 
area. The increased area is approximately 1.16 times the test section area. 

In this region test section flow and plenum flow merge to form diffuser inlet 
flow. The ejector action of the diffuser maintains the required plenum flov, 
through the flaps. During supersonic flow a shock wave is established in this 
region. 

For simulation purposes the above transition is assumed to occur at a line. 

The effective cross sectional area at this line being the numerical value of Ar 
used in the simulation. Equation (13) equates the momentum balance across this 
line. As in shown in this report equation 13 as well as the value of are the 
main significant factors in the characterization of the test section flow phenomena. 
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